Journal of Nonlinear Mathematical Physics 1998, V.5, N 3, 234-244. Letter 
I 

Similarity Reductions for 

a Nonlinear Diffusion Equation 

M.L. GANDARIAS, P. VENERO and J. RAMIREZ 

Departamento de Matematicas, Universidad de Cadiz, 
P.O. Box 40, 1510 Puerto Real, Cadiz, Spain 
E-mail: mlgand@merlin.uca.es 

Received November 19, 1997, Revised May 27, 1998, Accepted July 8, 1998 

Abstract 

Similarity reductions and new exact solutions are obtained for a nonlinear diffusion 
equation. These are obtained by using the classical symmetry group and reducing 
the partial differential equation to various ordinary differential equations. For the 
equations so obtained, first integrals are deduced which consequently give rise to 
explicit solutions. Potential symmetries, which are realized as local symmetries of a 
related auxiliary system, are obtained. For some special nonlinearities new symmetry 
reductions and exact solutions are derived by using the nonclassical method. 

1 Introduction 

Equation 



with n S Q\{0, 1} and A G M, corresponds to nonlinear diffusion with convection. When 
C = equation (|l]) becomes the one-dimensional porous medium equation 

«t = («")„ . (2) 

A complete group classification for (g) was derived by Ovsiannikov 23] and by Bluman 
E , Q . A classification for Lie-Backlund symmetries was obtained by Bluman and Kumei 
||]. The main known exact solutions of nonlinear diffusion (^) are summarized by Hill 
|15| . In [15, 16, 17], Hill et al deduced a number of first integrals for stretching similarity 
solutions of the nonlinear diffusion equation, and of general high-order nonlinear evolu- 
tion equations, by two different integration procedures. King [jn| obtained approximate 
solutions to the porous medium equation (g). 
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The basic idea of any similarity solution is to assume a functional form of the solution 
which enables a PDE to be reduced to an ODE. The majority of known exact solutions 
of (||) turn out to be similarity solutions, even though originally they might have been 
derived, say by a separation of variable technique or as traveling wave solutions. For 

u t = (« B )„ + C («»),, 

which is the Boussinesq equation of hydrology involved in various fields of petroleum 
technology and ground water hydrology, several exact solutions have been obtained by 
using isovector method |jj . 

More often than not the spatial dependent factors are assumed to be constant, al- 
though there is no fundamental reason to assume so. Actually, allowing for their spatial 
dependence enables one to incorporate additional factors into the study which may play 
an important role. For instance, in a porous medium this may account for intrinsic fac- 
tors, like medium contamination with another material, or in plasma, this may express 
the impact that solid impurities arising from the walls have on the enhancement of the 
radiation channel. Knowing the importance of the effect of space-dependent parts on the 
overall dynamics of the nonlinear diffusion equation, a group classification for 

«t = + /(*X«* + g(x)u m (3) 



was derived in by studying those spatial forms which admit the classical symmetry 
group. Both the symmetry group and the spatial dependence was found through consistent 
application of the Lie-group formalism. The behaviour of the interface of a related problem 
with (||) has been studied by Okrasinski [22]. 



In this work we use the invariance of equation ([l]) under one-parameter group of trans- 
formations to reduce the PDE (jl|) to various ordinary differential equations. Most of the 
required theory and description of the method can be found in g, |l|, |l|, ||, || . Following 
Hill |l7| , we have deduced some exact solutions of equation (|l[) by fully integrating 

the ODE's derived. 

An obvious limitation of group-theoretic methods based in local symmetries, in their 
utility for particular PDE's, is that many of these equations do not have local symmetries. 
It turns out that PDE's can admit nonlocal symmetries whose infinitesimal generators 
depend on integrals of the dependent variables in some specific manner. 

In H |(| Bluman introduced a method to find a new class of symmetries for a PDE. By 
writing a given PDE, denoted by R{x,t,u}, in a conserved form, a related system denoted 
by S{x,t,u,v} with potentials as additional dependent variables, is obtained. Any Lie 
group of point transformations admitted by S{x,t,u,v} induces a symmetry for R{x,t,u}; 
when at least one of the generators of the group depends explicitly on the potential; 
then the corresponding symmetry is neither a point nor a Lie-Backlund symmetry. These 
symmetries of R{x,t,u} are called potential symmetries. 

The nature of potential symmetries allows one to extend the uses of point symmetries to 
such nonlocal symmetries. In particular: Invariant solutions of S{x,t,u,v} yield solutions 
of R{x,t,u} which are not invariant solutions for any local symmetry admitted by R{x,t,u}. 
Potential symmetries for equation (0), when it can be written in a conserved form, have 



been recently derived [12] 
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In order to find potential symmetries of ([l]), we write this equation in the conserved 
form 

D X F - D t G. (4) 

The associated auxiliary system S{x,t,u,v} is then given by 

v x = (x + X)u, v t = (x + A) (u n ) x + (C - l)u n . (5) 

Suppose S{x,t,u,v} admits a local Lie group of transformations with infinitesimal generator 

Q Q Q Q 

x s = p(x, t, u > v )q^ + q( x , *i u > v )fo + r ( x > *) u , v )q^ + s ( x , *> u > v }fo- ( 6 ) 

This group maps any solution of S{x,t,u,v} to another solution of S{x,t,u,v} and hence 
induces a mapping of any solution of R{x,t,u} to another solution of R{x,t,u}. Thus 
defines a symmetry group of R{x,t,u}. If 

then @ yields a nonlocal symmetry of R{x,t,u}. Such nonlocal symmetry is called a 
potential symmetry of R{x,t,u} [||, ||. 

Motivated by the fact that symmetry reductions for many PDE's are known that are 
not obtained by using the classical Lie group method, there have been several generaliza- 
tions of the classical Lie group method for symmetry reductions. Bluman and Cole [3] 
introduced the nonclassical method to study the symmetry reductions of the heat equation 
and Clarkson and Mansfield |8[ presented an algorithm for calculating the determining 
equations associated with the nonclassical method. The basic idea of the method is that 
the PDE (||) is augmented with the invariance surface condition 

pu x + qu t - r = 0, (8) 

which is associated with the vector field 

d d d 

V = p(x, t,u)— + q(x, t,u)—+ r(x, t , u )^- ( 9 ) 

By requiring that both (||) and (|8|) are invariant under the transformation with infinite- 
simal generator (^) one obtains an overdetermined nonlinear system of equations for the 
infinitesimals p(x,t,u), q(x,t,u), r(x,t,u). The number of determining equations arising 
in the nonclassical method is smaller than for the classical method, consequently the 
set of solutions is in general, larger than for the classical method as in this method one 
requires only the subset of solutions of (|l]) and (|8|) to be invariant under the infinitesimal 
generator However, the associated vector fields do not form a vector space. These 
methods were generalized and called conditional symmetries by Fushchych and Nikitin || 



and also by Olver and Rosenau |2(j, 21] to include "weak symmetries", "side conditions" 
or "differential constraints". 

The nonclassical symmetries of the nonlinear diffusion equation (||) with an absorption 
term and f(x) = 0, as well as new exact solutions were derived in jL3|. In this work we 
obtain the special values of the parameter n such that nonclassical symmetries for (|l|) can 
be derived. We also report the reduction obtained as well as some new exact solutions. 
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2 Classical symmetries: Exact solutions 

In this case, we find that the most general Lie group of point transformations admitted 
by © is: 

3n + 1 

1. For n^O.l and C ^ we obtain a three-parameter group Q3. Associated with 

n + 1 

this Lie group is its Lie algebras £3, which can be respectively represented by the 
set of all the generators {l^}f =1 . These generators are: 

9 , , x d „ d Tr d u d /H „. 

_ V2 = {x + X) - + 2t - v 3 = -t- + — _. (10) 

3n + 1 

2. For n ^ —1,0,1 and C = , we obtain a four-parameter group Q^. Its 

n + 1 

infinitesimal generators {Vi}| =1 are: 

V( = Vi, v> = v 2 , v> = v 3 , 

{x + xf- c d 2(x + X) 1 - c d ( n ) 



vL 



-u- 



(C- 2)(C- 1) dx (C-l)(n-l) du" 
To ensure that an optimal set of reductions is obtained from the symmetries of (0), the 



optimal system is determined for A 7^ 0. The case A = was derived in [11]. In Table 1, we 
list the nontrivial optimal system {Ui\ with i = 1,2,3,4. We also list the corresponding 
similarity variables and similarity solutions. 

Table 1: Each row show the infinitesimal generators of the optimal system, the 
corresponding similarity variables and similarity solutions. 



% Ui Zi Ui 



1 Vi + aVi i.x + \) 2 a h(z){x + X)^ 

2 aV{ + V2 + 2V£ e- at {x + \) h(z)(x + A)^ 

2n 

h ^ eCl (x+X)^T 



3 aVi + VI ci (n -l)(x + Xf' 1 + log t 



(x + A)»+! 



tw a(C - 2)(x + A) 0-1 ..... . 

4 aV[ + Vl — /i(z)[(n - l){x + A)] 



In Table 1, case i = 1, the constant a G Q\ {2} and, in case i = 3, c\ 



n + 1 
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The ODE to which the PDE (@) is reduced by means of the infinitesimal generator U\ 

is 

dh ( 1 Ci+n + a-1 



dz\( a - 2) 2 h n ~ l n (« - 2) (n - 1) * 

ah{Ci-n + l) n-l fdh\ 2 d 2 h _ 

+ (a _ 2) 2 (n _ 1) v + -^v^J + ^"°' 

with 

C\ = Cn + an — C, 
after taking h = yn it becomes 



_dy / Ci + n + o - 1 _ 1 \ an (Ci - n + 1) y cPy _ 

^(a-2)(n-l)z (a - 2) 2 n J (a - 2) 2 (n - l) 2 z* dz* ' 

In particular, the second order nonlinear differential equation for y{z) obtained for i = 1, 
admits first integrals for some values of a. As an example, we consider 

a = 2 - 2n. 

The first integral for a = 2 — 2n is 

_y + Cy_ y_ | dy | V" _ 

z 2nz 2nz dz 4n 2 

where ci is the integration constant. If this constant is zero we obtain for a = 2 — 2n 

n 

(4n 2 +k 2 k 2 \ n-l 

(1 — n) z 4« 2 + (4kin 3 + kik2n — kik2) z^ I 

V = I P ' 

(4n 3 + k2n — k2) z 4 ™^- 4 ' 1 

with 

A: 2 = (2c- 2)n-4n 2 , 

(4n 2 +k 2 k 2 \ n-l 

(1 — n) z 4 « 2 + (4kin 3 + kik2n — k^) I 

" = -J- , k2 ' 

(4n 3 + k2n — k2) z 4 ^- 4 " 2 
Substituting in the similarity solution we obtain the exact solution 

(4n 2 +k 2 k 2 \ ' 

(n — 1) z 4 « 2 + (— 4kin 3 — kik2n + kik2) I 



i 

k 2 \ n-l 



(— 4n 3 — k2n + k2) z* 



When C = this is the well-known dipole solution 15]. 
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In particular, the second order nonlinear differential equation for h(z) obtained for 
i = 2, is 

dh(Cn + An-C a \ 2h(Cn + n-C + l) n - 1 / dh\ 2 d 2 h 

+ tz, t + — — 77o—„ + — ^- \-r) +^7 = 0. 



dz \ (n-l)z h n - x nz) (n - if z 2 h \dz J dz 2 

For the special value 

n — I 

we integrated once to obtain 

,„ ,dh h n (3n-l) h n ah 

h T z + f u + — = ci. 

dz (n — l)n n n 

After taking 

i 

h = y™- 1 

it becomes 

2y a a dy 

1 v ~r = ci, 

z nz z dz 
where c\ is the integration constant. If this constant is zero we obtain 

c a (n — 1) 

y = — — o — ' 

z z 2n 

so that 

= 2at 



ce 2at a(n-l)\ n - 1 



2 



h= {— 2 2^) X ~' 

Substituting in the similarity solution we obtain the exact solution 

/ ce 2at a(n-l)\^ 

The second order nonlinear ODE obtained, after taking h(z) = y l l n , for i = 3 is 
(n + l) 2 y~~ l e~ z dy d 2 y 2n dy n 2 y 



where c\ 



4ci 2 (n-l) 2 n 3 cte dz 2 n - 1 dz {n - if 
a 



n + 1 

The second order nonlinear ODE obtained, after taking h(z) = y l l n , for i = 4, is 



ete 2 \ 22: J dz 



n — 3 Q 

(n- l)»+i(n + l) 3 
where c\ = — 



16an 3 



240 



M.L. Gandarias, P. Venero and J. Ramirez 



3 Potential Symmetries 

In order to find potential symmetries of (|l|), we write this equation in the conserved 
form (H), where 

G = (x + X)u, 

F=(x + X) (u n ) x + {C -l)u n . 

d d 

The associated auxiliary system is given by (M). Besides X\ = — and X2 = — we obtain 

ot ov 

X3, . . . , Xq, given in Table 2. 

Table 2: Each row show the infinitesimal generators, the corresponding similarity 
variables and similarity solutions. 



(a) p q r s n C 



*2iU 2/2 1) 

X% x + A 7^ 0, 1 arbitrary 

n — 1 n — 1 

X 4 (x + A) 2nt -2u / 0, 1 arbitrary 

1— n 2/c 2n 3n + 1 

AT 5 fc(x + A)^ — (x + A)"^+T ^-1,0,1 ^— - 

n + 1 n + 1 

X 6 2(x + A)w -2(a + A) 2 n 2 - 2uv v 2 -1 -- 



X3, X4 and X5 project onto point symmetries of (JlJ) ; while X§ induces potential sym- 
metries admitted by ([[]). 

Solving the characteristic equation, we obtain the similarity variable z = t and similarity 
solution v = y/x + XE(t). In this case, E(t) satisfies E'(t) = 0, so E = const and we obtain 
the trivial solution 

Cx 



CiVx + A, 



u 



2(x + A)§ 



We must note that, although the infinitesimal p depends explicitly on v, the similarity 
variable does not depend on v. 



4 Nonclassical symmetries 

To apply the nonclassical method to (jl|) we require (|]) and (|8|) to be invariant under 
the infinitesimal generator (|9|). In the case g / 0, without loss of generality, we may 
set q(x,t,u) = 1. The nonclassical method applied to (jl]) give rise to four determining 
equations for the infinitesimals. 
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K du 2 dudx x + A du / 

- (n - 1) n^tt 11 " 2 + (n - 1) nra"" 3 - 2p^- = 0, 
ait a« 

^ <9 2 r 9 2 p nC op — nC \ n+1 
H dudx n dx 2 x + Xdx (x + A) 2 7 

-2 (n - 1) n%u n + ( 2^r - 2p^- - f^U 2 + (n - 1) pra = 0, 
ox \ du ox Ot J 

, d 2 r nC dr\ „ ( dr dp \ , . 9 

- 1 " 5? + —XBi) « + ( » + 2 ^ r ) " " ( " " '> * °' 

nC 

where f(x) = — — -. Solutions of this system depend in a fundamental way on the values 
of n. By solving the determining equations we obtain 

p = p 2 (x,t)u n +pi(x,t). 
We can distinguish now the following: if n G { 0, — 1, — - 1 we recover the classical sym- 



metries, and if n ^ |o, —1, — - j , we obtain that 

r = a 2 u n+2 + a 3 u n+1 + aiu™ +1 + — ^-r + a 4 u 2 + riu, 

u n 1 

where 

CP2 2p 2 2 1 d?>2 _2pi£2 

n(x + \y Q2 ~ (n + l)(2n + l)' ° 3 ~ n ax ' ° 4 " n + 1 

and pi, j»2 5 ^lj arid r 2 are related by two conditions. After considering the special values 
for n for which new symmetries different from Lie classical symmetries can be obtained, 
we can now distinguish the following: 

• for n — we recover the classical symmetries, 

• for n = ^ it follows that p 2 = and pi, r±, and r 2 are related by the following 
conditions 

api api 
-Pin-4 P1 — -2— =0, 

r ! d 2 pi 2 ng | 2 nC< -0 

ax dx 2 x + X dx (x + A) 2 

at ax 
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dt TlT2 dx 12 dx 2 x + A dx 

d 2 ?2 _ 2 nC dr 2 + 2 = Q 
dx 2 x + A dx 

We do not solve the above equations in general, but consider some special solutions: 

2 



1. Choosing p\ = k, r± = and C = 2, then r2 
nonclassical ansatz 

z = x — kt, u = I 



(x + A) 5 



and we obtain the 



fc(x + A)y ' 

where k ^ and satisfies the following ODE 
fc" + 2fc/i/i' = 0; 

whose solutions are 

fc 4 tanh(fc 4 (z + fc 2 )) , 7 - 

= ; , if kk 3 > and fc 4 = \/kk^, 



(12) 



fc 

fc 4 tan(fc 4 (z + fo)) 



, if fcfc 3 < and fc 4 = \J—kk%. 



This leads to the exact solutions 

'fc 4 tanh(fc 4 (2: + fo)) 



u(x, t) 



u(x, t) 



fc 4 tan(fc 4 (z + k 2 )) 



k(x + C) 



k k(x + C) 

2. Choosing pi = pi(x), T\ = r 4 (x) and r 2 = r 2 (x) we obtain 



if fcfc 3 > 0, 



, if fcfc 3 < 0. 



r(x, t, u) = —Ap'u + 



hp>{ C P [ + C 



Pi (s + A)pi ' (x + A) 2 
where p\ must satisfy the following equations 
C 3p[ p'{ ci 



2(x + A) 2 Pl 2p\ 2p 1 p' l ' 
2C 



x + A 



[2p?pr+4pipkf-(pi) s ] 



cpi 



(x + A) 2 



[cpi P '; + 8p 1P '; + 2(p[) 2 } 



C (3C - 2)pfpi + ^^(C - 2)Cp? - 5PM" 



(x + A) 3V ^ (x + A) 4 

-lQpiPiPl" + 3Qpi(p?) 2 - io(K)M = 0. 



We observe that these condition are satisfied forpi = fci(x+A) 2 if C €. ^ -, —1, 3 >. 
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For C = 3, pi = k\{x + A) we recover classical symmetries. 
5 

For C = — , we obtain 

1/3 4A; l n 



pi = k\(x + A) ' , r 



3(x + a) 2 / 3 ' 
Hence a nonclassical ansatz is 

Z = ^l( x + A) 2 / 3 -kt, u = (x + X)- 4 / 3 h(z), 

where h(z) satisfies the following ODE 
-2hh" + {h'f + -4kkfh 3 / 2 ti = 0, 
whose solutions are 

h 

Akfk 



= — ~r~oT tan(/cx A;4 (z + k 3 )) if kk 2 > 0, k^ = \Jkk 2 and k\ < 



= — ^— tanh(/jifc4(2 + /C3)) if /c£;2 < 0, /C4 = ^J—kk 2 and A4 > 
4k(k 

• For C = — 1, we obtain 

ki _ 4(/cin + 2/u) 

Pl ~^TA' r ~ (2 + A) 2 ' 

Hence a nonclassical reduction is 

z = - t, u = + A) 2 - 2) 2 , 

fcl fcf 



5j?i - 2/ci 



where satisfies ( [12] ) with fc = k\. 
3. Choosing ^2 = ^2 = and p\ = pi(x), we obtain r\ = —ip[ and / 



2pi 

where pi satisfies 

pfpi + 2^1^ ) 2 - 2kxp\p' x - k 2 p\ = 0, 



which is a classical reduction [11|. 



5 Concluding Remarks 

In this paper we have used the invariance of (|lj) under group of transformations to re- 
duce (|) to ODEs. We desired to minimize the search for group-invariant solutions to that 
of finding non-equivalent branches of solutions, consequently we have constructed all the 
invariant solutions with respect to the one-dimensional optimal system of subalgebras, as 
well as all the ODEs to which (JlJ) is reduced. For the equations so obtained, first integrals 
have been deduced which give rise to explicit solutions. Potential symmetries as well as 
nonclassical symmetries were used to obtain new solutions of (Q). The new solutions are 
unobtainable by Lie classical symmetries. 
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